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tract. Long mvalength aiioaallM in the 
total angnetlc fiald aaaaurad tqr Kagaat ovar the 
Unlt^ States and adjacent areas are inverted to 
an equivalent layer crustal ■agnetizatlon dlstri- 
bution. The a^el is based on an equal area 
dipole grid at the Earth's surface. Model resolu- 
tion, defined as the closest dipole spacing giving 
a solution having physical slg^flcance, is about 
220 lea for Magsat data in the elevation range 
300-500 kn. The nagnetizatliMi contours correlate 
veil with large-scale tectonic provinces. 

Introduction 

Techniques developed for interpretation of Pogo 
satellite aagnetlc antmaly data (Mayhew et al, 
1980; Mayhew, in press) were ained at 1) aodellng 
the anoealy field Itself and separating it from 
core and external fields, and 2) deriving aodels 
of large-scale crustal oagnetlzation which have 
geologic significance. These techniques have now 
been Applied to data taken over the U.S. by 
Magsat, which was considerably lower chan Pogo, 
and initial results are given in this paper. The 
approach has been to apply the well-known equiva- 
lent source technique (e.g. Oaspney, 1%9), adap- 
ted for spherical-earth geoatetry. The eqxiivalent 
source technique is a slaple method of generating 
a synthetic nagnetic anomaly field which fits the 
data over its full elevatlmi range. Magnetic 
dipoles are arrayed in an equal area grid at the 
Earth's surface and ■acheaatlcally oriented in a 
fixed direction. In previous applications, the 
direction has been taken to be cl» same as that of 
the core field. A sec of dipole «»iencs Is deter- 
ained by a generalized least squares method such 
that the dipoles collectively give rise to a mag- 
netic aniMuly field best fitting chat observed. 
Although previous applications have used only the 
anomaly in the total field as input data, in prin- 
ciple vector c<mponent measurements can be used 
directly )ust as readily. For the work reported 
here we have used total field anomaly data as c(»- 
puted from the measured vector components. The 
set of dipole moments is divided by a surrounding 
volume given by the dipole spacing and an arbi- 
trary vertical thickness, so Chat the resulting 
solution is expressed as apparent magnetization 
variation in a constant -thickness layer (thus, an 


lay«r wdcl)* This It In cimcrMt 6o « 
Mkgnntic Bodni for tho crust which would hswo con- 
stant ■sgnstisstion but vsrlsbls thlcknoss; s 
rssllstlc gsologlcsl aodsl would in ^Mrsl 
Include wsristions In both psrsastsrs* Msgnstlss- 
tl<m psrsMtsr solutions and ths sssoclstsd fit of 
aynthstic and aasaurad fialds display a charactar- 
is tic bahavlor as a function of dipola spacing. 
As tte apaci^ dacraasaa (ntmbar of ^pola aourcaa 
incraasas) thm fit of tha field iaprovaa rapidly 
but with dininishing ratums beyond a critical 
spacing. At about this spacing tha stability of 
tha solution, aa aaaaurad simply by tha paraaater 
standard deviation, begins to datariorata rapid- 
ly. Dacraaaing tha spacing further causes tha 
■agnatlsation values to taka on Increasingly 
larger positive and nagativa values. Tha critical 
spacing can be daflnad rather closely. For spac- 
Ings close to. but mit lass than, this Unit, 
awgnatlsatioa variations appear to be physically 
■sanlngful. in aa auch as they are systaaatic and 
correspond closely to large-scale geologic provin- 
ces (Mayhaw at al. 1980; Mayhav. In press). For 
snallar spacings. solutions are not physically 
Bsaningful. although they do give a slightly 
bettor fit to the data. Exanplaa of such " trade- 
off" graphs for Pogo data are given In tha above 
rafarancaa. A trade-off graph for Magsat data 
over tha western United Statas and ^Jacant araas 
is given In Figure 1. Tha stability Halt for 
Pogo data (above 450 ka altitude) corresponds to a 
source spacing of about 300 ka; for Hsgsat data 
(w>stly above 300 ka) we infer a Halt of about 
220 ka. Tlw lower data leads to batter resolution 
aodals. as expected. Mayhaw (In press) describes 
a aathod of datarainlng a aagnetlsation distribu- 
tion for a dipola array with spacing half that of 
tha critical spacing (for tha present data. 110 
ka). and this method has bean applied In producing 
tha aagnetlsation distribution described below. 
Equivalent layer magnetisation models can be read- 
ily transform^ to more realistic aodals of aagne- 
tlsation variation In a negnetlc crust of variable 
thickness, provided soma Independent inforaatlon 
is available for calibration. 


lUsults 


A virtue of the equivalent source techniqiM is 
that the anoeaiy field can be represented over a 
surface of constant altitude and displayed in aap 
fom. Figure 2 shoiis the equivalent source total 
field anoealy at 320 lea altitude computed from a 
136 km ** spacing equivalent source array. Data 
and sources within about 32* of the magneflc pole 
were not used in order to avoid external tlcld 
effects in the auroral sone. 

Figure 3 is a aagnetisatl<m model derived from 
Magsat total field an<MMly data* Contours are In 
tenths of A/m, and represent an apparent magneti- 
zation contrast distribution within an equivalent 
layer 40 km thick* The prlnlpal features are the 
same as in a model derived from Fogo data (Mayhew, 
in press), but the present model shows consider- 
ably more detail* Many of the magnetization 
anommlles correlate with large-scale tectonic fea- 
tures* In the western United States, the pattern 
corresponds to large-scale heat flow provinces 
(e*g* Lachenbruch and Sass, 1977), and may reflect 
Curie isotherm undulations* In particular, the 
Basin and Range ami Rio Grande Rift are regions of 
hl^ heat flow, ami are delineated by magnetiza- 
tion lows (possibly indicating a shallow Curie 
isotherm), while the relatively lower heat flow 
provinces of the Sierra Nevada and Colorado 
Plateau correspond to magnetization highs* The 
boundary between the Appalachian-Ouachlta belt and 
the Precambrlan cratlon is marked by a strong gra- 
dient along its langth* The Wichita Uplift, 
flanked on the north by the Anadarko Basin, is 
marked by a magnetization gradient* It is inter- 
esting that in the western midcontinent the boun- 
dary between Mesozoic/Cenozoic cover and older 
rocks marks a distinct change in orientation of 
the magnetization anomalies* The Reweenawan rocks 
of the Lake Superior Syncline appear to be associ- 
ated with a local magnetisation positive* Attent- 
ion is called to the similarity between the con- 
tour pattern of Figure 3 and the rsagnetic map of 
the United States based on filtered aeromagnetic 
data given by Sexton et al (this issue)* 
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Figure Captions 

Fig. 1. Trada-off aa a function of dipole spacing 
(kilomatars) batwaan 1) the fit of tha aquivalant 
source aagnatic an<»aly field to the field obser- 
ved at satellite elevations (tenths of nT) and 2) 
"stability'* of inversion as indicated by standard 
devia(;lon of solution parameters (tenths of A/m). 
Optimal solution Inferred to be at source spacing 
of about 220 km. SO means standard deviation. 

Fig. 2. Equivalent source representatlcm of magne- 
tic anomaly' field at a height of 320 km. Source 
spacing is 136 km. Contour interval InT. Albers 
equal area projection. 

Fig. 3. Apparent magnatizatlon contrast in a 40 km 
thick layer. Distribution is obtained by inver- 
sion of Magsat total field anomaly data. Model is 
sarlaa of staggered 222 km - spacing dipole grids 
computed separately, but machine contoured togeth- 
er without smoothing on 111 km grid. Contour 
intarval is 0.1 A/m. Dashed lines indicate gener- 
alised tectonic province boundaries which are geo- 
graphically related to tha magnetization distribu- 
tion. Letters represent tha following; S: Sierra 

Nevada block, Nt: Basin and Range province, CPt 

Colorado Plateau, R: Rio Grande Rift, SR: Snake 

River Plain, P: waatam boundary of Great Plains, 

A: Appalachian Ouachita front, W: Wichita up- 
lift, R: Kentucky anomaly, LS: Lake Superior 

Synclina. Boundary batwaan Masosoic/Canozoic 
cover and older rocks la indicated by short dashed 
line. Albers equal area projection. 
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